Bethe Ansatz solution of a new class of Hubbard-type models 
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We define one-dimensional particles with generalized exchange statistics. The exact solution of a 
Hubbard-type Hamiltonian constructed with such particles is achieved using the Coordinate Bethe 
Ansatz. The chosen deformation of the statistics is equivalent to the presence of a magnetic field 
produced by the particles themselves, which is present also in a "free gas" of these particles. 



Studies on strongly correlated one-dimensional (ID) 
systems have deeply influenced modern concepts in 
many-particle physics. These systems may be of strictly 
linear extension like quantum wires, quantum-Hall bars, 
(quasi-) ID organic metals and spin chains, or higher di- 
mensional systems where rotational degrees of freedom 
are not relevant, as in the Kondo problem Exactly 
solvable models play a crucial role in the theoretical 
description of ID systems, since there, most of the re- 
sults obtained by approximate methods are not reliable 
(especially at low temperature) . They also serve as a reli- 
ability test for various approximations, which are applied 
to models in higher dimensions. 

An important example is the Bethe Ansatz (BA) solv- 
able lD-Hubbard model (HM) j|. The HM, originally 
proposed in Ref . [Q , describes electrons hopping on a D- 
dimensional lattice, while experiencing an on-site interac- 
tion. The two-dimensional HM is believed to capture im- 
portant features of high-T c superconductivity || . Away 
from half filling, the 1D-HM exhibits Luttinger Liquid 
(LL) behavior ||. 

Various generalizations of the HM have been proposed; 
most of them emerge from changing explicitly the form 
of the Hamiltonian 0). The coupling of the fermionic 
degrees of freedom to a phononic bath has been studied 
in Ref. Recently, Schulz and Shastry || proposed a 
class of solvable Hubbard-type models in which the hop- 
ping amplitude of the particles with a given spin orienta- 
tion is modulated by a gauge field (entering the Hamil- 
tonian via a Peierls-like substitution) , which depends on 
the density of particles with opposite spin orientation. 
A similar gauge-like potential has been employed to cou- 
ple several Hubbard chains (h]] , where "M-colored" elec- 
trons in M chains perform intra-chain hopping, which is 
(gauge-like) affected by the density of electrons in the 
other chains p0[ . The effective inter-chain interaction 
resembles an electromagnetic field due to the motion of 
electrons moving along each chain. The models studied 
in Ref. || show LL behavior even in the absence of the 
interaction: the asymptotic behavior of correlation func- 
tions depends on the coupling strength between the gauge 
field and the particles, which is responsible for the change 
in periodicity of the wave function. The gauge potential 
breaks the time reversal symmetry of the Hamiltonian. 
Charge and spin excitations have anyon-likc character, 



and the scaling properties deviate from the "fermionic" 
LL theory g. 

Another way of extending a model, is to modify the "na- 
ture" of the particles entering the Hamiltonian, without 
changing its formal structure jnj. In the present paper 
we choose this approach to extend the HM: we keep the 
formal structure of the Hubbard Hamiltonian unaltered, 
but we change the particles from fermions to particles 
whose wave functions obey a generalized exchange sym- 
metry. We call the statistics of such particles Deformed 
Exchange Statistics (PES) E|- Here, a variant of the 
DES introduced in Ref. JT3J is used. The corresponding 
particles constitute a non-Abelian realization of the sym- 
metric group S N (see Ref. fjjj). The "deformed" ID HM 
is solved exactly by means of the Coordinate BA (CBA) 
and the Bethe Equations (BE) are compared with the 
results in Refs. § and @. 

We define the HM Hamiltonian for particles obeying 
DES as 

H = -t (ftji+i,* + h-c) + U £ v lA v lA , (1) 

i,a i 

where f i;Cn f} a (a G {|, j}, respectively a G {1/2, -1/2}) 
obey the deformed relations 

f},ijfk,a' + Qj'k fk,cr'fj !(7 = $j,k Sera', ^ 
fj,erfk,tr' Qk, 'j fk.o' fj,o — 0, 

Qff = {Ql'f)- 1 = QK'f , (3) 

[fl,Q?£] = \h«,Q°%'] = o (4) 

Relations (|^) are formally analog to quon commutation 
rules ]T^| , but here the deformation parameter depends 
on indices (j, a\k, a'). The operators Vj jCT = fjafj,c 
are the particle-number operators. Equation (||) guar- 
antees that the particles are representations of Sn, 
whereas Eq. (|]) ensures the standard commutation re- 
lations [vi,<T,Vj,<T'] = 0, [vi,<r,fj t<T '] = Si,j5a,cr>fl a ,, and 
[ v i,mfj,a'] — —Si.jfia.cr'fj.o-', providing a well defined Fock 
representation of the algebra defined in Eqs. (||). We 
point out that DES is defined for any Qj'^ fulfilling (||), 
(^); we choose the deformation parameter depending on 
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spins and coordinates, but also such that the solvability 
of the model (|]) is ensured: 



-j,k 



o i-0(cr+o-')gi4>-(cr— a') 
„i*-(<T— cr') 

— iip(a+a' ) p i& ■ (a— a') 



j > k, 
3 = k, 
j < k. 



(5) 



For fixed j, k, cr, a', is a C-number. 

The deformation defined above becomes integrable 
if <f> - $ ± ip(0) = mod 2n holds. This implies 
($ - $;V(0)) = (0;0) V (vr; tt) mod 2?r. We note that 
inserting definition (0) in Eq. (0) implies the Pauli prin- 
ciple since Qj'j 7 = 1, and that two particles with different 
spin on the same site obey deformed anticommutation re- 
lations since Q]j ± 1. Equation (|) is trivially fulfilled 
(q is a C-number) as well as the consistency relation (||) 
together with Eq. (||) . The fermionic case is obtained set- 
ting $ = $ = 0; ip = 0, whereas the hard-core bosonic 
case is covered via $ = $ = 0; ip(±l) = n. Furthermore 
we mention that for spin- 1/2 particles and factorizing 
spin- and coordinate dependence of the deformation, ev- 
ery spin dependence in the deformation can be written 
in the linear form used here. 

It is worth mentioning that (||) implies fixing an order 
on the ring. Such an order distinguishes one (given) lat- 
tice site from the others in that it defines the "beginning" 
of the ring, where one starts counting. Periodicity means 
not only that the wave function has to be periodic, but 
simultaneously demands the result being independent of 
where the starting point was set. What is equivalent, is 
fixing a period Pq = { j\ , . . . , ji + L } on the infinite 
periodic chain. Consistency of the Periodic Boundary 
Condition (see below) with this induced order is given if 
the results are independent of Pq . In the sequel it will be 
seen that this condition is fulfilled. 

Now we show that the Hamiltonian ([!]) is exactly solv- 
able by means of CBA. We apply Periodic Boundary 
Conditions (PBC), which means fj+L,a = fj.a- The gen- 
eral iV-particle state in a chain with L sites can be written 
as 



I*) = 



E *Oik)^Oik)|o>, 



(6) 



where we have used the abbreviations j 7 = 

(j 7 (l),- ; -,j 7 (JV)), ^Hh k) = ■■//,(„>, ^(m,; 

for arbitrary 7 G Sn', 1 is the identity in Sn- For in- 
distinguishable particles in D ^ 2, a Fock base state is 
uniquely determined by the coordinates and the spin- 
configuration; it must not depend on the order in which 
the particles are created. This requires 

*C?7 It ^)^{h It t) = *Cv IV ° ^^Uy W ° 7r ) 

(7) 

for any 7, 7' S Sn, since changing the order in creating 
particles means permuting both the coordinate and spin 



indices of the operators // , „ , in ■ Equation (0) 
fixes the symmetry of the wave function ^ . In the case of 
fermionic statistics, ^(j-y I7071") = sign(7)\&(ji | vr) , where 
sign(7) = +/— for even/odd permutations. Every per- 
mutation can be decomposed into a product of neighbor 
exchanges Xk = k «-> (k+ 1). Taking 7 = Xn ■ ■ ■ °X2 °Xi 
as a simple and relevant example to fix the idea and fur- 
ther defining the partial products as 7/ = Xi° ■ ■ - °X2°Xi 
for I < n, we obtain 



*(i 7 17° 



n 

-)=n[ 



-Q 



x'(l),x(l) 



k), 



(8) 



1=1 



where x{l) = j 7l07 r(i), s(l) = cr 7i07r(/) , x'(l) = j 7l07 r(;+i), 
and a' (I) = cr 7l0w(i+1 y. 

Due to the DES (0), special care must be taken when 
employing the symmetry of to restrict the coordinate 
space to a sector of ordered coordinates (called "order re- 
gion" in the following) . An order region can be character- 
ized by a permutation t £ Sn which connects the actual 
set of coordinates (ji, . . . , jn) with the ordered one so 



that (ji, ...,j N ) = (h 



(i)' 



r(AT) 



), where i\ < 



< ir. 



The corresponding wave function, which is defined in this 
order region, is labeled by r: *ff T (jn' \n) = ^i(j TO n' |ro7r). 
At first, it is necessary keeping this label for writing 
the Schrodinger Equation, since the hopping term can 
connect different order regions. Due to the Pauli prin- 
ciple, only two order regions are needed to solve the 
Schrodinger Equation (let these two order regions be con- 
nected by Xk)- For 7 = Xk, Eq. (0) simplifies to 



k) = -q: 



l^r(i k), 



(9) 



where x = v (fc ), a = a nh) and x' = V(fe+i), cr' = ov (Jk+1) . 
In the special case of a doubly occupied site, iu = ik+i, 
the symmetry relation on the border of an order region 
leads to 

*x* (?k) = \Xk o tt) = -Ql'^tCi k), (10) 

Equations (^) and ( |l0| ) allow writing the Schrodinger 
Equation in terms of a single ^ T [16|. Since each ^> T 
can be obtained from using Eq. (ph, we consider ^1 
in the following. The label 11 will be omitted, and the 
coordinate vector j is chosen being ordered. We insert 
the Bcthe wave function 



*(ik) = E ^(^k)exp 



N 



JmPiT'(m) 



into H\^) = E\*ff). For pairwise distinct coordinates 
we obtain the energy E in terms of the momenta pf. 
E = — 2t 53;=i cospi. Its form is unaltered by the de- 
formation parameter q. Using Eqs. @ and ([To|), the 
Schrodinger Equation corresponding to doubly occupied 
sites explicitly contains the deformation parameter, and 
the scattering matrix S reads 



2 



S(X,X';m) 



i{X - A')e i *-^-"')P m , m+1 



u 

21 



i(X - A') 



j7 
21 



(ID 



where we have defined A = svn{p n /f m -\), A' = 
sin(p 7r / (m+1) ), = ^(m), and // = cr ;r(m+1) . The per- 
mutation operator P m , m +i is defined by its action on 
the amplitudes Pm,m+i A(ir'\n) = A(ir'\xm ° tt)- 

The scattering matrix S 1 fulfills the Yang-Baxter Equa- 
tion (YBE) 

S(A, A'; m + 1) S(A, A"; m) 5(A', A"; m + 1) = (12) 
S(X', A"; m) S(X, A"; m + 1) S(X, A'; m), 

where A" = sin(p 7r '( m+ 2))- F° r the validity of the YBE, 
it is important to notice that P and ((J, — //) do not com- 
mute, since fi and // contain the spin permutation, which 
is affected by P: P mi „ l+ i(^-/i') = (jt/-^)P m)m+1 . Ap- 
plying PBC, which also demands independence from the 
pre-chosen order, we obtain 

• •, JjvK) = QpBC^Uk+i, ~,3n,3i +L, ..J k +L\%), 

(13) 



where tt is the cyclic permutation tt 
the function Qpbc reads 



/ 1 ... N-l N \ i 

U+i...fe-iJi and 



Qpbc = (-l)^" 1 ) J] J Q 



m=l 1=1 
k 

= (-1)W-Vq% >0 l[Q B ^ lmV (14) 

m— 1 

where Qi?^ = ^[(^-l)^)^*^] is the spin depen- 
dent part, and Qb.o = 1 the spin independent part of 
Qpbc @|- Thus, the deformation affects the BE in two 
distinct ways: (a) via the prefactor Qp B c by imposing 
PBC; its form arises from the deformed symmetry of the 
wave function induced by the algebra (||); (b) due to a 
modified S'-matrix in case of on-site deformation. The 
BE are 



tt ijsinpj - Ca) - Tt 

QbaQb.o i{smpj - Ca) + Tt ' 



Qs l 



(15) 



N 

n 



T *(Ca " C&) + § 



A Z ( Sin P' - Ca) - Tt 



n 



= l *(Ca - Cf>) - § <2^T i(sinp; - Ca) + Tt ' 

where the factor Qs = e -4 * arises from the S matrix in 
diagonalizing of the transfer matrix of the "spin prob- 
lem", Eq. (|l3|). Using the definitions of Qb,o,Qb,o, Qs, 
and the integrability condition $ — <t±ip(0) = mod 2ir, 
the factors due to the DES are 



Qs l 



= e -i(JVT-l)^(l) 



Qb,iQs 
Qba 



= e t[(JVi-i)tK-i)-(JV T -iW>(i)] < 



(16) 
(17) 



We see how the above contributions involve only the spa- 
tial part of the deformation parameter; a purely spin- 
dependent deformation of the statistics does not affect 
the spectrum of the model. Even in this case, how- 
ever, the eigenstates are different due to the modifica- 
tion in the S matrix. This can also be made trans- 
parent by representing the /, y by fermionic operators 



n- 



In a forthcoming pa- 



per |18J we will describe in more detail the connection 
between fermionization and DES, as well as the connec- 
tions to Refs. |) 14 1 . Here we mention only that, defining 
<t>cr = —(N a — l)-0(2er), the phases appearing in the BE's 
decompose into <£>| and <Fr — $|, found in Ref. [Q. In 
contrast to the result in [||, the factors N a — 1 appear. 
The reason for this is, that the phases in Ref. |j| were 
produced by particles having opposite spin, whereas here 
the relevant phase comes from particles having equal spin 
projections. Thus the particle feeling that phase has to 
be excluded. 

From the BE's (|l5|), the total momentum can be ex- 
tracted: 



N 



, V>(2<r) 
L 



(18) 



From this it is seen that the total momentum, though 
periodicity is implied, is not necessarily a multiple of 
2tt/L, which is the case for the undeformed model. This 
is not surprising, since the spectrum is that of a twisted 
fermionic model. The total momentum can be obtained 
in shifting iVj momenta of the undeformed model, which 
become 2im/L for U — > 0, by — (iVj — l)tp(l) and the re- 
maining N± momenta by — (iVj — l)ip(— 1). That meant 



N„ 



2tt. 



(19) 



with $ CT as defined above in comparison with Ref. Ill ] 
and li a , for fixed u, being distinct integers modulo L. 
Expression ([l9|) shows the physical meaning of DES: The 
deformation of the particles' statistics is equivalent to a 
magnetic field generated from the particles themselves. 
Such a magnetic field depends on how many particles 
are in the system. In particular (|l9| ) shows how the en- 
ergy of a gas of free (U — > 0) particles is not simply the 
sum of the single particles' energies, but it describes a 
system of interacting particles. Such interaction purely 
comes from the deformed statistics. However, the non- 
interacting limit is to be taken carefully here. This can 
be understood by facing the second BE. In order to keep 
the phase, two or more spin rapidities are forced to cither 
coincide with each other or with one or several sin(pj). 
Such contributions have to be linear in U /At. 
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In conclusion, we have extended the Bethe ansatz for 
particles with deformed exchange statistics or, which is 
equivalent, for solutions with generalized exchange sym- 
metry. A method already introduced in Ref. |[3| is here 
generalized for problems including inner degrees of free- 
dom like spin. Using this technique, a new class of 
generalized Hubbard models could be shown being ex- 
actly solvable. The spectrum equals that of a fermionic 
model with spin dependent boundary phases, which in 
fact turns out being a general feature of integrable de- 
formed exchange statistics. The eigenfunctions however 
are different due to the non-fcrmionic exchange symme- 
try. A preliminary study showed that these phases in 
general already contribute to the thermodynamic limit. 
We will discuss the details elsewhere fis}] . 
A systematic study of deformed models is most interest- 
ing for at least two reasons. First, "integrable deformed 
statistics" can be an alternative way of handling com- 
plicated interactions of integer statistics' particles. In a 
forthcoming paper Jl8| we will review these results from 
another view-point; namely fermionic models with corre- 
lated hopping. There, a general statement on integrabil- 
ity of deformed exchange statistics will be drawn as well 
as for models with correlated hopping. 
Second, it might be relevant for modelling the edge states 
in fractional quantum Hall effect physics, the chiral Lut- 
tinger liquid behavior of which has recently been ques- 
tioned 
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